Observable form of pulses emitted from relativistic collapsing objects 
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In this work, we discuss observable characteristics of the radiation emitted from a surface of a 
collapsing object. We study a simplified model in which a radiation of massless particles has a 
sharp in time profile and it happens at the surface at the same moment of comoving time. Since the 
radiating surface has finite size the observed radiation will occur during some finite time. Its redshift 
and bending angle are affected by the strong gravitational field. We obtain a simple expression for 
the observed flux of the radiation as a function of time. To find an explicit expression for the flux we 
develop an analytical approximation for the bending angle and time delay for null rays emitted by 
a collapsing surface. In the case of the bending angle this approximation is an improved version of 
the earlier proposed Beloborodov-Leahy-approximation. For rays emitted at i? > 2Rg the accuracy 
of the proposed improved approximations for the bending angle and time delay is of order (or less) 
than 2-3%. By using this approximation we obtain an approximate analytical expression for the 
observed flux and study its properties. 



PACS numbers; 04.40.Dg, 95.85.-e, 97.60.-s 



I. INTRODUCTION 

Neutron stars and black holes are compact relativistic 
objects. Light propagating in the vicinity of such ob- 
jects is affected by a strong gravitational field. For the 
description of the photon propagation under these con- 
ditions the general relativity is required. This is evident 
for black holes where the gravitational potential at the 
horizon (in the units) is of order of 1. Though a value 
of the gravitational potential at a surface {Rns ^ SM©) 
of neutron stars (M^s ^ 1.5Mq) is smaller and is of 
order of ~ 0.2, it was demonstrated recently that general 
relativistic effects might be important for understanding 
the features of the radiation coming from these objects 
0-0. In particular, one of the explanations of X-ray 
bursts is connected with possible explosions of the mat- 
ter which is accumulated on the surface of a neutron star 
entering in a binary system as a result of "sucking" the 
matter from a usual-star companion. It was shown (see 
0, 0, 0) that the profile of the X -ray-curve observed 
by a distant observer can be explained only if the effects 
of general relativity are taken into account. In partic- 
ular, according to the general relativity, because of the 
gravitational bending of light rays, a distant observer can 
see a part of the opposite side of the neutron star which 
is invisible in a fiat spacetime. The radiation emitted 
from this part gives an important contribution and has 
a visible impact on the form of the signal from the X-ray 
burst. 

The effects of the general relativity also modify consid- 
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erably light curves for continuous in time radiation from 
the surface of a collapsing star as seen by a distant ob- 
server. These effects were studied in detail by Ames and 
Thorne 7], Jaffe and Lake 0. In this study main 
attention was focused on the details connected with light 
emitted near (unstable) circular photon orbits at 3rg/2, 
where Vg is a gravitational radius of a collapsing star. In 
such considerations there were usually adopted a number 
of simplifying assumptions, such as: (1) Spherical geom- 
etry; (2) Dust-like (pressure free) equation of state; (3) 
Radiation comes only from the (free-falling) surface of 
the star; and (4) It is continuous in time. 

In this paper we would like to discuss a slightly differ- 
ent set up, when the assumption (4) is violated. Namely 
we assume that the radiation emitted from the surface 
of a collapsing spherically symmetric stellar object has a 
profile of a sharp in time pulse. Such radiation may occur 
in different situations. For example, suppose a neutron 
star or a proto-neutron star looses its stability as a result 
of the accretion of matter onto it or due to the softening 
of equation of state at the center which is supposed 
to be already several times higher than the nuclear den- 
sity . During the collapse, the matter density of a com- 
pact object is growing and the whole system evolves into 
the much higher density than the normal nuclear den- 
sity beyond which new hadronic phase transi- 
tions might take place One might expect a possible 
sharp-in-time emission of massless particles (photons and 
neutrino) during such phase transition 

Since at the moment there are a lot of theoretical un- 
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certainties in the estimations and detailed description of 
these processes, we shall study a simplified (toy) model. 
Namely we assume that a radiation of massless parti- 
cles has a sharp in time pulse profile and it happens at 
the surface at the same instant of time (from a point 
of view of a comoving observer). The time required for 
the radiation to reach a distant observer depends on the 
position of a radiative region on the collapsing surface. 
For this reason the pulse emission results in a continuous 
flux received by the observer during some finite interval 
of time. During this interval the flux as well as the red- 
shift factor changes. In principle, knowing the redshift 
and light curves allows one to obtain direct information 
about the collapsing body at the moment when the radi- 
ation occurs. We shall demonstrate that if the emission 
occurs at the moment when the radius of a collapsing 
body is i? > 4.5Af, the observed characteristics of the 
pulse radiation with a very high accuracy can be written 
in an explicit form in terms of elementary functions. We 
shall use this result to study how different (gravitational 
redshift, Doppler and focusing) effects affect the form of 
redshift and light curves. 

In section 2 we consider null rays emitted from a sur- 
face of a spherical collapsing object. The equations of 
motion of the freely collapsing surface and for rays prop- 
agating in the surrounding Schwarzschild geometry are 
well known. In order to fix the notations we collect the 
required formulas in section 2 in the form they are used 
in the further sections. The time of arrival of null rays 
emitted by a collapsing surface and their bending angle 
are given by elliptic integrals. In section 3, we propose 
an explicit analytic expressions which approximate these 
quantities with very high accuracy. For rays emitted at 
R > 2Rg this accuracy is of order (or less) than 2-3%. 
This approximation is used to obtain the expressions for 
the intensity and flux for a short flash observed by a dis- 
tant observer (section 4). The results of numerical cal- 
culations of the redshift, bending angle and light curves 
are given in section 5. In section 6, the results are sum- 
marized and discussed. 



the dust-like equation of state so that a solution for the 
interior of the collapsing object is a well known Tolman 
solution while the surface follows a radial geodesic 
in the Schwarzchild geometry [l^ . 

Denote by r the proper time as measured by an ob- 
server comoving with the collapsing surface. We suppose 
that the collapse starts at r = and the initial surface 
radius is i?o- Then, in the parametric form (using the 
conformal time parameter 77), the equation of motion of 
the surface is 



R 



Rq 



(1 -I- cos 77), 



T = \ — 2. (77 + sin 77) 



(3) 
(4) 



In the same parametric form the Schwarzschild time i*^*^^ 
corresponding to the proper time r at the surface is 



2Af In 



yjRo/2M - 1 + tan(77/2) 



v/i?o/2A-f - 1 - tan(77/2) 



+2A/Vi?o/2Af - 1[77 + ^(77 + sin 77)] . (5) 



We assumed that &'>{t = 0) = 0. 
We denote 



(6) 



This quantity has the meaning of the invariant radial 
velocity. Then in a general case the four- velocity of the 
collapsing surface is 



where 



v^' = {d&^/dT,dR/dT,0,0), 
d&'> 1 



(7) 



(8) 



II. BASIC EQUATIONS 

A. Freely-falling spherical surface 

We consider a photon emitted from a collapsing spher- 
ical surface and propagating to the observer at inflnity 
in the background of Schwarzschild metric 



ds^ = -fdt^ + f-^dr^ + r'^dn , 
/-/(70-l-2M/r, 



(1) 
(2) 



where M is the mass of the collapsing object. We adopt 
the natural units, c — G — h ^ 1. InQ and ^ the 
motion of a spherical surface during the gravitational col- 
lapse was discussed under assumption that the dynam- 
ical role of the pressure can be neglected. We also use 



— = Vif{R)—— = ^ 



(9) 



Using eq.|(21) and lO, we obtain the following expres- 
sion for the velocity Vi of a freely falling surface with the 
initial radius i?o at the moment when its radius is R 



2M y/1 - R/Ro 



We also have 



R ^1 - 2M/Ro 



dt^^'> _ Vl - 2M/Ra 



dr 1 - 2Af/i? 



dR 
d7 



-^2M/R^l- R/R 



■ 



(10) 



(11) 



(12) 
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Propagation of null rays emitted from the 
surface 



Consider a photon emitted from the surface. Its tra- 
jectory hes in the plane. Without loss of generality we 
assume that it coincides with a plane of the fixed coor- 
dinate (/), so that the vector of the 4-momentum of the 
photon is 



(13) 



Because of the symmetry of the Schwarzschild metric, 
E = —pt (the energy at infinity) and L = pg (the angu- 
lar momentum) are constants of motion. Instead of the 
angular momentum L we shall use the impact parameter 



I = L/E. 
The radial momentum p^ is 



P 



aEZ . 



where 



Z = Z{1, r) = y/l-l^f{r)/r^ 



(14) 



(15) 



(16) 



Here and later a denotes a sign function which takes the 
values -|- and — for a forward (pT > 0) and backward 
{p"^ < 0) emission, respectively. 

There is an upper limit for the impact parameter, Imax, 
for a photon which being emitted from a surface of radius 
R reaches the infinity: 



R 



VfW) 



(17) 



Suppose a null ray emitted by a collapsing surface reaches 
the infinity at 9 = and has the impact parameter 
Then the emission angle, /3, as measured by an observer 
comoving with the surface, is defined by the relation 



cos Pa {I, R) 



(7Zil,R) - 
1 -aviZ{l,R) 



(18) 



We consider only the light emitted with /? < 7r/2. A 
tangentially emitted photon with respect to a comoving 
observer(/? = 7r/2), for a collapsing surface [vi < 0), is 
possible only for a backward emission. The correspond- 
ing impact parameter I — It is determined by the condi- 
tion 



Z = -Ui . 
Solving this equation one obtains 

R 



h 



Vl - 2M/Ro 



(19) 



(20) 



For R = 3-\/3\/l — 2M/RqM the parameter It is equal 
to 3-^/3, the value corresponding to the (unstable) circular 
photon orbit. To escape delicacies connected with more 



complicated behavior of the photons near this orbit we 
assume that R > 3-^/3^/1 — 2M / RqM . In this case for a 
photon which reaches the infinity the possible ranges of 
an impact parameter are < Z < Imax and Zt < ^ < Imax 
for a forward and backward emission, respectively. A 
discussion of the allowed ranges of the impact parameter 
for the smaller radius up to i? ~ 2M, can be found in 
and 1. 



Redshift 



Let p^if^ be 4-momentum of a photon emitted by a col- 
lapsing surface which has 4- velocity t)^'^)''. Then v'^'^^ = 



-Ph V 



is the frequency of the photon as measured by 



a comoving to the surface observer. Since an observer 
at infinity is at rest, its 4-velocity is — Sq and the 

observed frequency of the photon is i^*-"^ = — 
where is the 4-momentum of the photon at infinity. 
For a given ray, the redshift factor $ is defined as the 
ratio of the emitted frequency v^''^ to the observed at 
infinity frequency v^"^ 



(21) 



For a ray with the impact parameter I emitted from the 
surface of the radius R one has 



l-aViZ{l,R) 



For a freely falling surface we get 



y/l - 2M/Ro 

1 - 2M/R 



(22) 



(23) 



'(2M/i?)(l - R/Ra 
1 - 2M/Ro 



■Z{l,R) 



D. Bending angle 

We use the freedom in the choice of spherical coordi- 
nates to put the angle 6 in the direction to an observer 
at infinity to be equal to zero, 0^°^ = 0. Consider a null 
ray emitted by the collapsing surface when its radius is 
R and which reaches the distant observer. Suppose its 
impact parameter is I. Then such a ray is emitted by the 
collapsing surface from the region at the angle 6^*^^ . For 
forward emission this bending angle is 



3(e) 



= e(;,i?). 



Q{l,R) = I 

Jr. 



dr 



r^Z{l,r) ■ 



(24) 



(25) 
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For a backward-emission a photon before it reaches the 
infinity should pass through a turning point, rt < R, 
which is determined by Z{1, rt) = 0, or 



r2/(l-2M/rO - f 



(26) 



One can see that, ioi I — I 
as expected. Then we get 



= i?V(l-2M/i?), n^R 
2Q{l,rt)-Q{l,R). (27) 



III. ANALYTIC APPROXIMATION 

A. Approximation for bending angle 

It is convenient to use the dimensionless quantities 

x^M/r, q = M/R, i=l/M. (34) 
We also denote 



E. Arrival time 

Consider a null ray with the impact parameter I emit- 
ted from the collapsing surface at the moment t when 
it has the radius R{t). Denote by t^^^ the time when 
it reaches a distant observer at radius for the for- 
ward/backward ray. One has 

dr 



so that 



Z = Z{l,q) = dl-Pq^l-2q), 



Zil, R) = z{i, M/R) . 



(35) 



(36) 



Using these notations the integral 8, (|25|l . can be writ- 
ten in the form 



f{r)Z{l,ry 



(28) 
(29) 



e 



I 



dx- 
Z{l,x) 



(37) 



- Ml 



dr 



dr 



As already was mentioned, for a given surface radius R 
there exists a critical value ^max 



f{r)Z{l,r) Jnir) f{r)Z{l,r) 



It is evident that t*-"' oo when rg — > cxd. For this rea- 
son it is more convenient to consider a finite quantity, 
the time difference between arrival of two null rays emit- 
ted at two different moments of proper time, r and Tg, 
respectively. For the second ray, emitted at Tg, we put 
I = 0. Such a ray goes radially. We denote this time dif- 
ference by A<(/;r, Te). In the limit when r*-"^ — s- co this 
quantity remain finite. For the forward ray it is given by 
the following expression 

r, Te) = t^'^ (r) - (Te) + T(/, i?(r)) 



i/(gyT 

(30) rays emitted from this surface must have 
a fixed q and I G [0, Z„ 
to 0. 

In a flat metric with / 
integral (pTzIl analytically to get 



- 2q) and the 

^ ^max ■ For 

the function Z changes from 1 
1, one can calculate the 



e 



flat 



arccos( Y 1 — Pq'^). 



(38) 



i?(re) -i?(r) + 2M In 



RjT,) - 2M 
R{t) - 2M ' 



Leahy % discovered that in a wide range of its arguments 
the exact integral for the bending angle can be approx- 
imated by a simple analytical expression. A simple ele- 
gant form of the approximative expression was proposed 
later by Beloborodov . In the notations adopted in the 
present paper this formula reads 



where 



T{l,R) 



dr 



Z{l,r) 



Similarly for the backward ray one has 
At_(/; T, Te) = t^^\T) ~ &\t,) + 2r(/, rt) 



(31) 



(32) 



1 - cos /3 = (1 - cos6i('^))(l - 2MIR) 



(39) 



where (3 is defined by H18() . It is easy to show that this 
relation corresponds to the following approximation for 




r(/,i?(rXJB3) 



9n = arccos 



Z -2q 
l~2q 



(40) 



+R{t^) + R{t) - 2rt + 2M In 



{R{t)-2M){R{t,)-2M) 
{rt - 2MY 

The integrals for O and T (see relations (|25|) and (|32|l . 
respectively) can be expressed in terms of the elliptic 
functions. However, for practical calculations it is very 
convenient to have approximations for these quantities in 
terms of simple elementary functions. In the next section, 
we develop high accuracy analytic approximations, for 
the integrals Q[l,R) and T{l,R). 



We shall refer to this relation as to Beloborodov-Leahy 
(or BL-) approximation. 

The typical accuracy of the BL-approximation is of 
order of 1% for the light rays emitted from the surface 
R = 6M. For this reason it works well in the study of 
light emitted from a surface of a neutron star In a 
recent paper Q BL-approximation was used for fitting 
the light curves of X-ray pulsars. For smaller R the ac- 
curacy of the BL-approximation is worse. For example 
it becomes of order of 10% for R = AM. In order to 



5 



use this approximation for our purposes we first slightly 
modify it to improve the accuracy. 

To improve the BL-approximation, let us compare ex- 
pansions of 8 and 80 at small values of I. We have 

e = g/ + ^(2~3g)q3f3 



— (21-70g + 60(72)g5[5 



(41) 



3 ;'3 



Qo = ql + j^i^ - 3q)q' I 



5 t5 



(. . . denotes terms of higher order in L) This gives 

We write our ansatz as follows 

'z{L,q)-2q 



Q{i,q) = arccos 



where 



l-2q 



l~Zil,q) 



For small I 
Z 

We write 65 as follows 



l-2q 



\q^P + Lx2-iq)qH^ + . 
I lb 



224 



(42) 



(43) 



(44) 



(45) 



(46) 



(47) 



For (3 = 1 this expression correctly reproduces the expan- 
sion of 8(/, q) for small q up to the terms of the order of 
0{q'^). Numerical calculations show that the accuracy of 
the approximation is very good for the following choice 



/3 = 3.5, 55 = 0.0884. 



(48) 



Thus the approximate expression for the forward-emitted 
rays is 



^eil,q = M/R), 



(49) 



where 9 is given by (|^ . l|17|l . and The relative 

error Ag = (6'jf'' — Q)/6^'^^ of the approximate expression 
is very small. It is less than 0.5% for R > 4.5Af for all 
the allowed values of I. Figure shows the quantity D = 
10^ X Ae as a function of s = l/lmax{q), where q = M/R. 
The parameter s is chosen so that for every value q and 
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^Ae (Ag is a relative error 
. The plots are shown for 



FIG. 1: The figure shows D = 
for 0) as a function of a — I /I 
the following 4 values of R/AI — 4.5, 5, 6, 7 (the smaller value 
R/M the higher is a plot in the vicinity of s = 1). The curves 
are practically indistinguishable from D = line for s < 0.5. 
For this reason we plotted only s > 0.5 parts of the graphs. 



the impact parameter I from to its maximum allowed 
value Imax the value of s belongs to the same interval 
[0,1]. For R = AM the error Ag is slightly larger. It 
is still less than 0.8 % every where excluding a narrow 
vicinity of s = 1 where it reaches 2%. 

We shall use the formula (|44|l to approximate the bend- 
ing angle for the forward emission. For the backward 
emission the approximate formula is 



g(e) 



[l, R) « 29(f, M/n) - 9(f, M/R) , 



(50) 



where rt is defined by H2t)|) . 



B. Approximation for arrival time 

Now we consider the arrival time. Using the dimen- 
sionless version of T, T = T/M, we can rewrite the ex- 
pression l|32|) in the following form: 



T = Til,q) 



Z{i,x){l + Z{i,x))' 



We want to obtain an analytic approximation for T. 

Let us first assume that the function / which enters 
Z changes slowly, and put / =const in the integral H51|l. 
This integral can be calculated exactly 



To = 



1 - \li-rwj_ 

fq 



(52) 
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We restore the dependence / on q and use this expression 
with / = 1 — 2(7 as a starting point for our approximation. 
The corresponding expression can be written as 



q 



(53) 



One can check that this approximation is very good for 
smah q. 




FIG. 2: The figure shows D — 10^ At (At is a relative error 
for T) as a function of s = l/lmax- The plots are shown for 
the following 4 values of R/M — 4.5, 5, 6, 7 (the smaller value 
R/M the lower is a plot in the vicinity of s = 1). The curves 
are practically indistinguishable from D = line for s < 0.5. 
For this reason we plotted only s > 0.5 parts of the graphs. 



To improve it, let us compare expansions of T and To 
at small values of I. We have 



i<Zp + i(g3-2g4)P 



2 lb 



This gives 



1 

16 



q't + . . 



(54) 
(55) 

(56) 



Our ansatz for the improved approximation is 

f = To + Q(Z) . (57) 

Using 146(1 one can conclude that in order to provide 
the correct q'^r behavior of the correction at small I the 
function Q must have the following expansion 



Q{Z) = 



We write our ansatz as follows 
1 



Q{Z) = ^Z^ + a^Z^ + aiZ^ . 



(59) 



Numerical calculations show that the accuracy of the ap- 
proximation is very good for the following choice of the 
parameters 

03== 1/15, 04 == 1/25. (60) 
Thus the approximation for the time delay takes the form 



r=- + Q(z), g(z) = iz2 + lz3 + lz^(6l) 

q 4 15 25 



til, R) = Mf{I, M/R) 



(62) 



The relative error At = {T — T) /T of the approximate 
expression is again very small. It is less than 0.5% for 
R > 4.5M for all the allowed values of /. Figure |21 shows 
the quantity D — 10'^ x At as a function of s = l/lmax- 
For R — 4Af the error At is slightly larger. It is still 
less than 1 % every where excluding a narrow vicinity of 
s — 1 where it reaches 3%. 

For the calculations of the intensity of the light from 
a radiating body observed at infinity we shall need the 
expression for the derivatives Tg and Tj. The first one 
can be easily found from the definition of T l|51|l 



Z{l,x){l + Zil,x)) 



(63) 



For the derivative we shall use the approximate ex- 
pression obtained by the differentiation of T given by 
(|61|l with respect to I. The approximation for the deriva- 
tive works slightly worse than the approximation for 
T. We denote 



(64) 



the relative error. The figure 13 shows D — 10^ At, as a 
function oi s — 1/ Imax ■ The maximum value of the error 
(for R — 4.5M) is near s — 1 and it reached 5%. 



IV. A SHORT FLASH FROM A COLLAPSING 
SURFACE 

A. Flux and intensity for a short flash 



We shall use superscripts (e) and (o) for emitted and 
observed radiation, respectively. To characterize such a 
radiation one uses the notions of the flux and intensity. 
The fiux, F, is a power of the radiation per unit area 
emitted by a radiating surface or received by the detec- 
tor. The intensity, /, is a power of either emitted or ob- 

(58) 

served radiation per unit area per unit solid angle. Quite 
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FIG. 3: The figure shows D = 10^ At, (At, is a relative 
error for Tj) as a function of s = l/lmax- The plots are shown 
for the following 4 values of R/M = 4.5, 5, 6, 7. In the region 
where D < 0, the smaller value R/M the lower is a plot in the 
vicinity of s = 1. The curves are practically indistinguishable 
from D = line for s < 0.5. For this reason we plotted only 
s > 0.5 parts of the graphs. 



often the spectral characteristics, called the specific flux, 
Ty, and specific intensity, 2^, are also used. They are 
"per unit frequency" versions of F and /, respectively. 
Evidently, one has 



F 



dvTi, , I 



dvl,. 



(65) 



The intensity / and specific intensity depend on the 
angle 9 between the direction of radiation and the normal 
to emitting (observing) area. The flux and specific flux 
are obtained by integration of these quantities over the 
directions 



F ■ 



dm , 



dilXij 



(66) 



where dfl is an element of a solid angle. For an isotropic 
case with azimuthal symmetry. 



dn = 2tt sin OdO . 



(67) 



The observed intensity 7*^°^ depends on the direction, 
or (in our case) on the angle 9 between the normal to the 
detector surface and the direction of observation. Instead 
of the direction angle 9 one can use the impact parameter 
I of an incoming photon. For large rg one has 9 = l/ro, 
and d^^°^ = 2T:ldl/rQ. Thus the relation can be 
written as ^3 



Idl- 



.27r/('') 



(68) 



In the integrand, the intensity I^"-* is now considered as 
a function of impact parameter I: 1^°'^ = 

The impact parameter I is conserved along the photon 
trajectory. It specifies the trajectory itself. Thus for a 
given value / one can determine both, the bending angle 
6'('^^ (the angle on the surface of a collapsing body where 
the ray was emitted) and the angle /3 at which it was 
emitted. 

One can show (see e.g. Exercise 22.17 in ^3)' ^^8.t 
the quantity I^{v)/v^ remains constant along a photon's 
world line. This quantity is proportional to the num- 
ber density in phase space for photons and its conserva- 
tion follows from the Liouville's theorem (see e.g., and 

E^)- . . . . 

Consider a light ray with the impact parameter / emit- 
ted from the collapsing surface at the moment of the 
proper time r, and let t be the time when it reaches a 
distant observer at tq. For fixed tq relations (|28|l and 
(|30|l can be used to determine r — T{t, I). 

The specific flux as measured by a distant observer at 
time t is 



(69) 



Here $ = $(l, i?) is given by ea. (|22|l . and R is the radius 
of the collapsing surface at the moment when a null ray 
with the impact parameter / leaves it in order to reach the 
observer at tq at the moment t. We denote by t the cor- 
responding proper time on the surface. For a fixed time 
of arrival, t, and the moment of emission r, one can use 
the relations H28() and (|30() to determine the correspond- 
ing value of the impact parameter / = /((r). Thus the 
integral over I in eg. 1)69(1 can be rewritten as an integral 
over the proper time, r: 



. ^ 27r /• 



dTW<i>'^ljfHl,iye,T), 



where 



W = l 



dl_ 
d7 



(70) 



(71) 



The expression for W can be obtained by differentiat- 
ing eq. and eq. (|^ with respect to r for a fixed value 

of i^"'* . Using the analytic approximations for the arrival 
time, we get (for details see Appendix) 



W = J-^ \Z-avAR- 



Here dt^'^^ /dr is given by ((11(1 . and 
J+ = l + ^Q\Z) , 



J-=J+ 



AMRZ {rt + MQ'^J 
Tt ~ 3M 



(72) 



(73) 



(74) 
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Here Q',^ = Q'if-Hn)), Z = Z{l,R), Z = and 
<i> — ^{l,R), and the arguments in these functions are 
I = lt{T) and R = R{t). We also use a notation Q' = 
dQ/dZ. A transition from the forward to the backward 
ray regime occurs at the point where Z = that is at 
I = Imax)- It is easy to see that at this point W- = W+, 
as it is expected. 

We shall need an explicit form of W for a special case 
when a ray is emitted radially and hence 1 = 0. In this 
case only forward emission is possible and we have J7+ = 
1 . Thus the corresponding value Wq takes the form 



Wn = 



2M y/1 - R/Ro 



RyJ\ - 2M/R 
1 - 2M/Ro 



(75) 



It is also easy to check that for a tangentially emitted ray 
(I = It) W vanishes. 

For a very short in time flash from the surface at the 
moment Tg, the intensity can be approximated as 



(76) 



In this case the integration over r in (|70|l is trivial and 

we get for the observed specific fiux the following 
expression 



2tt 



(77) 



Here R^ = i?(Te), 1^ = /((rg), = W{le,Re): and 
<I)g = $(;e,i?g). The factors We and <i>e, which enter 
this relation, do not depend on the frequency. Thus for 
the observed flux we have 

F(°)(t) = / dv,Ti°\t) = I dve^-'Til\t) . (78) 



(79) 



Using relations 1()5|) and 177|l one obtains 

F(°)(t)^^Ty,$-4/(e)(;^)^ 

where the intensity is 



(80) 



It should be emphasized that the simplicity of the ex- 
pression l|7i^ is a consequence of the assumption that 
the pulse has a (5-like form. In a general case the in- 
tegration over the time and frequencies cannot be per- 
formed explicitly. For this reason, a so-called monochro- 
matic approximation is often used in the discussions of 
the continuous radiation from relativistic objects (see e.g. 
0). We note also that in the adopted approximation all 
the expressions in H79|) are defined in an explicit form in 
terms of elementary functions. This radically simplifies 
the study of the light curves. 



B. Normalized flux 

We denote the flux registered by the distant observer 
at the moment t as F^°\t) and we denote by F^°\0) 
the flux at the moment when the first ray arrives to the 
distant observer. Such a ray is emitted forwardly (a = + 
and Z = 0) at = 0. It is convenient to normalize 
the observed time-dependent flux F^°^{t) to the value 
F(°)(0). We denote this ratio 



F(°)(i) 



F(°)(0) ■ 

Using the relation H79|) one obtains 



Wo 



$0 



/(^)(0) 



(81) 



(82) 



(83) 



The explicit expressions for the 'first-ray' quantity Wq is 
given by (|75|l . The equation H23(l determines $e and <I>o. 
In the latter case one needs to put a = + and Z = I. 
Relations (O, GSJ), and allow one to find We- 

Let us discuss now the last factor which enters the 
expression for F(<). In general, the intensity of the radi- 
ation from the surface of star can be written [23 as 



=a + bcos{l3{l)), 



(84) 



where a and b depend on the details of the emission pro- 
cess. In this work, we calculate two extreme cases: (A) 
a ^ 0, 6 = 0, and (B) a = 0, b ^ 0. The former cor- 
responds to the isotropic emission from a optically thick 
object and the latter corresponds to Lambert's law. For 
I = one has (3{l) = and l'^'^' = a + 6, so that 



a + 6cos(/3(/)) 
a + b 



(85) 



In the first case Ia = 1, while in the second one Ib = 
cos(/3(/)). 



V. REDSHIFT AND LIGHT CURVES FOR A 
FLASH FROM A COLLAPSING SURFACE 

To illustrate the obtained results, we consider now spe- 
cial examples. As a first example, we consider a neutron 
star which looses its stability. In this case an initial ra- 
dius Rq is Rns = 12 — 20 km and the mass is of order 
of M - 1.5Mq [23, and hence Rq/M = 5.4 - 9. An- 
other example is a proto neutron star RpNS ^ 20 km 
and M ~ 1.5Mq [jj]. In this case Rq/M = 9. In this 
section, we discuss in detail two cases with initial radii: 
Rq = 5AM and Rq = 9M. 

For a freely falling surface with Rq = 5.4A/, the turn- 
ing point rt on the trajectory of a backward ray lies 
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within the vahd range of the analytic approximation, 
rt > 4.5M, provided Re > 4.8Af. In accordance with 
this we choose Re = 4.8M (case I). 

For a freely falling surface with Rq = 9M, the ana- 
lytic approximation can be applied to the emission at 
Re > 5.5M. In this case we calculate a bending angle, 
redshift and a fluxes registered by a distant observer for 
the following 3 values of Re/M — 5.5, 6.5, 7.5 (cases Ila, 
lib, and lie, respectively). 



A. Arrival time 

Since the arrival time depends on the position of the 
radiating region on the surface, a light emitted at the 
same moment Te from the surface Re reaches a distant 
observer through some flnite interval of time t. Dur- 
ing this period, the observed brightness and frequency 
are changing. First ray arriving to the distant observer 
is emitted from 0e — and it has I — and a — +. 
The time when this first ray reaches a distant observer 
depends on the position vq of this observer. This time 
becomes infinite when tq — s- oo. It can be made finite 
if, for example, instead of t one uses the retarded time 
u = t — r* , where r* is the tortoise coordinate. We choose 
u — for the moment of arrival of the first ray. At given 
radius Am = At. The maximum arrival time difference 
is assigned for the backward ray emitted with an impact 
parameter It, and it is 



Atmax — A.t-{lT]Te,Te) 

+2Re - 2rt + 4M In 



2T{lT,rt)~T{lT,Re), 
{Re - 2M) 



(86) 



in ~ 2M) 

In the case I, for Rq/M = 5.4 and Re/M = 4.8, the 
time delay is calculated is Atmax/M — 13.8. In the case 
II, for Rq/M = 9 the time delay for different values of 
Re is given in the Table. 

Time delay for Ro /M = 9 



Case 


Ila 


lib 


He 


R,/M 


5.5 


6.5 


7.5 


Atmax/M 


16.9 


15.4 


14.4 



In what follows it is convenient to use a normalized 
arrival time difference defined as 6 = At/Atmax- We 
shall call this quantity the time parameter. The time 
parameter is always changes in the interval [0,1]. The 
time parameter for forwardly emitted light increases as 
I increases from ^ = to Imax- The backward emission 
starts with Imax and ends at It and the time parameter 
for a backward emission is increasing as I changes from 



B. Bending angle and redshift 



is shown in Fig. ^ (an upper curve) . For a comparison a 
similar function for a static surface (of the radius Re — 
4.8M) is shown at the same Figure (a lower curve). 




FIG. 4: Bending angle as a function of the time parameter 
S. The upper curve shows the bending angle for a freely- 
collapsing surface for Rq = 5.4M and Re — 4.8A/. The lower 
curve show the bending angle for a static surface with Re = 
4.8M. 



Figure |S1 demonstrates the bending angle as a function 
of the time parameter for 3 different cases Ila, lib, and 
lie. The smaller is Re/M, the faster is the radial motion 
of the radiating surface, and the larger is the observed 
region with the backward emission. As a result the range 
of bending angle for smaller values of Re/M becomes 
larger, as it can be seen at the Figure 




FIG. 5: Bending angle as a function of the time parameter S. 
The three graphs show the bending angle for the same initial 
radius _Ro — 9M for 3 different values of R^: (a) 5.5M, (b) 
6.5M, and (c) 7.5M. 



The bending angle as a function of the time parameter 
is a monotonously increasing function. For the case I it 



The frequency observed at infinity is different from the 
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5 5 

FIG. 6: Redshift factor for a freely collapsing surface as a FIG. 7: Light curves for 1 = 1 (lA) and I = cos/3(/) (IB) 
function of the time parameter 6 for the cases I, and IIa,b,c. for a freely-collapsing surface with Re/M = 4.8 and Ro/M = 

5.4M. 



frequency at emission because of two reasons: (1) Differ- 
ence of the gravitational potential at the point of emission 
and observation, (gravitational redshift), and (2) The ve- 
locity of the emitting surface (Doppler shift). The pho- 
tons emitted from the surface of experience the same 
gravitational redshift independent of their angular posi- 
tions (bending angle) of emission. However Doppler shift 
depends on the relative velocity of the surface of emis- 
sion with respect to the distant observer, and hence it 
depends on the bending angle (or the impact parameter 
I). The corresponding total frequency shift can be calcu- 
lated from ea. (|23|l . Since the arrival time depends on the 
impact parameter as well, the frequency shift then can 
be plotted as a function of the arrival time. The calcu- 
lated ratio of emitted frequency to the observed one, $, 
for a short flash as a function of the time parameter S is 
shown in Fig. |S1 Three curves which meet one another 
at (5 = 1 correspond to the three cases IIa,b,c. The forth 
curve corresponds to the case I. 

It is interesting that the redshift due to the gravity 
is substantially cancelled by the Doppler shift for the 
tangentially emitted light. Really, using (fT^ and we 
obtain for the redshift factor $ the following expression 



$7 



l-vf 



V fiRe) 

For a free fall from the radius Rq one has 



(87) 



o 1 - 2M/R 
l-vf ^ 



Thus 



$7 



1 - 2M/Ro 



^1 - 2M/Ro 



(89) 



It means that for the "last rays" (that is for rays with 
I — It), the redshift depends only on the initial radius 



Rq and does not depend on the radius of emission R^. 
For this reason the three curves IIa,b,c in Fig. |B| merge 
at the same value 3/V7 ~ 1.134 at (5 = 1 (that is for 
I = It)- Relation (|89|l also shows that for Rq = oo the 
gravitational redshift is exactly cancelled by the Doppler 
shift i]. 

Since for direct radial rays {I = 0) both effects "work" 
in the same direction, one can expect that for a given Rq 
the redshift will be larger for smaller values of R^- The 
Fig. IHl clearly demonstrates this. 



C. Light curves 

Let us discuss now normalized flux as a function of the 
time parameter 6. We call the corresponding graph a 
light curve. In case I {Rq = 5M, Re = 4.8M), the light 
curves for a short flash from the collapsing spherical sur- 
face are shown Fig.Q. The plot lA is a light curve for 
1=1 and the plot IB is a light curve for I = cos (3(1) 
(Lambert's law). Similar light curves for a static spheri- 
cal surface are shown for comparison at Fig. |S1 

For the case II {Rq — 9M) the light curves for A and 
B type of the radiating surface are shown in Fig. |51 and 
Fig. 1101 respectively. Each of the figures contains 3 curves 
corresponding to IIa,b,c cases. 

Let us discuss now qualitative behavior of the light 
curves. The observed normalized flux F{t) given by eq. 
(|82|l is a product of 3 factors: (1) a kinematic term 
Weit)/Wo, (2) a redshift factor ($e(^)/*o)~^ and (3) 
a normalized intensity of the emission I. The third fac- 
tor depends on the model of the radiating surface and it 
does not depend on the arrival time. The first two fac- 
tors are time dependent. The arrival time dependence of 
We is essentially determined by the factor of |Z — avi\ 
in ea. H72l) . which is a decreasing function of S and van- 
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FIG. 8: Light curves for I = 1 (lA) and I = cos/3(0 (IB) for 
a static surface with Re/M = 4.8. 




5 

FIG. 9: Light curves for 1 = 1 for the cases IIa,b,c. 



ishes for (5 = 1. Hence every light curves should cross the 
zero-flux axis at i5 = 1. For a static surface Vi — and 
Z (and hence We{t)) vanishes at (5 = 1, where I = Imax- 
For a collapsing surface Vi < Q and the observable flux 
vanishes not for Imax (where Z = 0) but for the backward 
emission with I = It- Hence one can expect longer dura- 
tion of observed flux for the emission from a collapsing 
surface compared to the emission from a static surface. 
The effect of motion of the collapsing surface becomes 
stronger for larger Vi. For example for a given i?o = 9A/, 
^tmax is calculated to be larger for smaller for which 
Vi is larger (see Table). 

The redshift factor $ depends basically on the relative 
receding velocity of the emitting region (determined by 
the bending angle) with respect to the distant observer. 
The relative receding velocity is decreasing as the bend- 
ing angle is increasing. Since the arrival time difference 




FIG. 10: Light curves for I = cos /3(Z) for the cases IIa,b,c. 



becomes larger for a ray with larger bending angle, one 
can expect the enhancement of a factor, for larger 

5. The main effect of the frequency shift for the observed 
flux due to the collapsing surface is the enhancement of 
the flux for lately arriving rays. As a result the shape 
of the light curve for a collapsing surface is changing 
from that of a static surface in such way that the flux 
decreasing in 5 is delayed and the sharp forward peak at 
6 — becomes a rather smooth peak as shown in Fig. jT)). 
For sufficiently large collapsing velocity, for example for 
Rq = 9M and Re = 5.5M, one can observe that position 
of the peak in the light curve also changes from (5 = 
to a later arrival time 5^0 for the isotropic intensity 
profile (A) as shown at Fig. Q. The emission angle with 
respect to the normal to the surface, /3, varies from 
to 7r/2 as the impact parameter I varies from to Imax 
and further to It- Hence the intensity profile of (B) with 
= 6cos(/3(/)) suppresses the enhancement due to the 
factor <I>~^ for lately arriving rays substantially as shown 
in Fig. ^1 The light curves from the static surface for 
intensity profile (A) and (B) are also shown in Fig. |Slfor 
comparison. 

VI. DISCUSSION 

In this work, we discussed characteristics of the radi- 
ation emitted from a surface of a collapsing object. We 
studied a simplified (toy) model in which a radiation of 
massless particles has a sharp in time pulse profile and it 
happens at the surface at the same instant of time (from 
a point of view of a comoving observer). In this approxi- 
mation both integrals over the time and frequency which 
enter a general expression for the observed flux can be 
taken. As a result we obtained an expression for the nor- 
malized flux as a function of the time of the observation. 
We demonstrated that for a short in time flash, the ob- 
served normalized flux can be expressed as a product of 
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three terms: a kinetic term (W), a redshift factor ($), 
and the intensity of the emitted radiation. The intensity 
of the emitted radiation is model dependent. In par- 
ticular it depends on the model of the radiating surface 
region. The other two factors are universal and depend 
only on the equation of motion of the collapsing surface. 
The dependence of W and $ on the arrival time is quite 
different. We assume that the collapse starts at some ra- 
dius Rq and the emission occurs at the radius Re- Under 
this assumption all the characteristics of the radiation 
such as its duration, redshift, and the form of the light 
curves (for chosen intensity of emission) depend only on 
these 2 parameters. To obtain the light curves one needs 
to integrate the equations for the light propagation in the 
Schwarzschild metric. The corresponding integration can 
be performed in terms of elliptic integrals. But the calcu- 
lation of the flux requires solving an inexplicit equation 
which contains the elliptic integrals. 

In order to solve this problem we developed an analyt- 
ical approximation for the bending angle and time delay 
for null rays emitted by a collapsing surface. In the case 
of the bending angle this approximation is an improved 
version of the earlier proposed BL-approximation [^|^. 
For rays emitted at i? > 2Rg the accuracy of the approx- 
imation for the bending angle and time delay proposed 
in the present paper is of order (or less) than 2-3%. By 
using this approximation we obtained an explicit formula 
for the observed normalized flux, which not only allows 
one very efficient numerical calculations of the character- 
istics of the radiation, but also appears to be useful for 
understanding the qualitative features of the radiation. 

Even for a static object the effects of the General Rel- 
ativity allows one to "see" a part of its opposite side 
surface. For a collapsing object this effect is more pro- 
found. As a result, the duration of the flux is elongated. 
Another difference is that the sharp decrease in time for 
the static surface is delayed and "smoothed out" so that 
the peak becomes broader. For a sufficiently large col- 
lapsing velocity the peak position can even be shifted to 
S > 0. We demonstrated these features by considering 
two examples of collapses starting at Rq = bAM and 
Rq = 9M. 

Though in this paper we focused on a model of brief 
in time flash emission, some of its results (improved an- 
alytic approximation) might be of the interest for other 
astrophysically interesting problems. 



APPENDIX A: CALCULATION OF dl/dr 

Let us fix Te in H31() and consider time delay A+(Z; r, t^) 
as a function of / and r. Consider a rays which reach the 
observer at the same time t. For these rays A+(Z;t, Te) 
is the same. This establish the relation between I and r: 
I = lt{T). To obtain this expression we differentiate H31f) 
with respect to r, keeping t fixed. 

For a forward ray one has 



fZ^ dl 



(Al) 



We denote by a dot the derivative d/dr. Using (O we 
get 



I 



(Z - Vi)V'^ 
Z{dT/dl) ■ 

To calculate the partial derivative 

dT _ dT dZ 

'dl ~ dz'm 

we use the approximate expression (|61l) . We have 
dT 

— = R + MQ\ 



dZ_ 

'dl 



R^Z ■ 

Combining these results we obtain 

R{Z - Vi)i(^^ 



l\ 



where 



J+ = l + iM/R)Q' . 



(A2) 



(A3) 



(A4) 



(A5) 



(A6) 



(A7) 



For backward rays the calculation of I is similar. Using 
the equation (|33|l for rays reaching the observer at the 
same time t one obtains the following equation for ^((r) 



tie) 



R 



f{R)Z(l,R) 
Using we find 



dT ■ 
'dl 



1 + 2 



Tin) 



rt 



fin) 



0(A8) 
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n = 

We also have 

zin) = 0, 



llPin) 
in - 3M) 



Zin) = i/fin). 



(A9) 



(AlO) 



Using the approximate expression for T, (|61|l . we have 
T., =Z(rO + [n+MQ;j^. (All) 



Here Q'^^ ~ Q'if ^(^t))- Using these relation we obtain 



Tin) 



= -II- 



2M 



fin) ^\l{n-iM) 

This aUow us to rewrite (|A8|) in the form 

AMRZ n + MQ[ 



[n + MQ'] . (A12) 



i^^^Z + v,)R = li 
Finally we obtain 

W- = II 



J+ 



n - 3M 



R{Z + Vi)f'' 



(A13) 



(A14) 



J- =J+ 



AMRZ n + mq;^ 

^ n - 3M 
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